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After reduction techniques, two-loop amplitudes in M = 4 super Yang-Mills theory 
can be written in a basis of integrals containing scalar double-box integrals with rational 
coefficients, though the complete basis is unknown. Generically, at two loops, the leading 
singular behavior of a scalar double box integral with seven propagators is captured by a 
hepta-cut. However, it turns out that a certain class of such integrals has an additional 
propagator-like singularity. One can then formally cut the new propagator to obtain an 
octa-cut which localizes the cut integral just as a quadruple cut does at one-loop. This 
immediately gives the coefficient of the scalar double box integral as a product of six tree- 
level amplitudes. We compute, as examples, several coefficients of the five- and six-gluon 
non-MHV two-loop amplitudes. We also discuss possible generalizations to higher loops. 
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1. Introduction 



Recently there has been renewed interest in the perturbation expansion of Af = 4 super 
Yang-Mills. This was motivated by the discovery of a twistor string theory |]J that captures 
the perturbation theory of the maximally supersymmetric Yang-Mills theory (pMSYM). 
Twistor string theory has opened new avenues and has inspired new ideas for the com- 
putation of tree level amplitudes of gluons [@SB^,[7| ! ^,E30E'E30E3'E^O=E1JI31 
and one-loop amplitudes of gluons in QCD 0J21|,|22|, TV = 1 |§,0,|^,|^|^,|^,|2| and 
Af = 4 [|30|,|31|.[3^,|33|,|341,p5| super Yang-Mills. Before twistor string theory was introduced, 



the study of pMSYM at one-loop was mainly motivated by two facts: one is the decom- 
position of a QCD amplitude, A^ CD , with only a gluon running in the loop in terms 
of supersymmetric amplitudes and an amplitude with only a scalar running in the loop, 
A scalar , (see []36| for a review), 



^qcd = A 



7V=4 



a aN=1 I /iscalar 
4/1 chiral + A 



Xl) 



where A^ =4 has the full Af = 4 multiplet in the loop and A^r^ only an Af = 1 chiral 
multiplet. The other motivation is a surprising proposal of Anastasiou, Bern, Dixon, and 
Kosower (ABDK) that two- (and, perhaps, higher-) loop amplitudes in pMSYM can be 
completely determined in terms of one- loop amplitudes ||37|| . This idea was inferred from 
studying collinear and IR singular behavior of the higher loop amplitudes. The conjecture 
is given in terms of normalized 2-loop amplitudes Mn = A n /A^ ee and in dimensional 
regularization, as follows 
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\ (MW(e)) 2 + /(e)M<?>( 2e ) - \u + O(e). 



;i.2) 



This relation was explicitly verified for four-gluon amplitudes in |37| (see also section 7 
of ]38[]). Also based on collinear limits ||39|| , the schematic form of a relation analogous 
to ( |1.2| ) was proposed for higher loops |4(J. Very recently, an explicit formula, analogous 
to (|1.2|), for the three- loop four-gluon amplitude was obtained and successfully verified 
in ||41|| . It is the aim of this paper to make some modest steps towards the calculation 
of higher loop amplitudes in pMSYM. The main motivation is to prepare the ground for 
future tests of the ABDK proposal. A proof of ( |1.2j ) would lead to the solution of pMSYM 
at two loops as a general solution to the one-loop problem can be obtained in terms of 
tree- level amplitudes by using quadruple cuts 0]. This is possible thanks to the cut 
constructibility of one-loop amplitudes in pMSYM proven in [[12] and the decomposition 



1 



in terms of scalar box integrals, with rational functions as coefficients, also given in E!| . 



See also [p0|,|3l|,|32],p3,35l for other techniques in pMSYM at one loop. 



At two loops, a similar decomposition in terms of some given set of integrals is ex- 



pected by using Passarino-Veltman or similar reduction procedures ||43|| . Unfortunately, 
the complete basis of two-loop integrals is currently unknown. However, scalar double 
box integrals are a natural ingredient of such a basisB. In this paper, we concentrate on 
the calculation of the coefficient of certain classes of planar scalar double box integrals. 
These are the integrals that arise in scalar field theory with a massless scalar running along 
internal lines and with the double-box structure depicted in fig. 1. 



(a) (b) 

Fig. 1: The two possible different structures of planar scalar double box integrals, 
(a) Double boxes. (6) Split double boxes. Note that the momenta of the external 
lines is given by the sum of the momenta of external gluons. 



The momenta of the external legs in fig. 1 are given by sums of momenta of external 
gluons. 

We propose a method for computing the coefficient of any scalar double box integral 
given in fig. la when at least one of the two boxes has two adjacent massless three-particle 
vertices. We also give the form of the coefficient of any double box given in fig. lb. In 
order to distinguish between the double boxes in fig. la and in fig. lb we refer to the 
former simply as "double boxes" and the latter as "split double boxes" . 

Our original motivation was the successful use of quadruple cuts in the calculation 
of one- loop Af = 4 amplitudes [j34| . The basic idea is that at one- loop only scalar boxes 



contribute [fl2fl . A quadruple cut singles out the contribution of a given scalar box and 
localizes the integration over the loop momentum. The combination of these two facts 



1 In fact, the four-gluon amplitude is given only in terms of scalar double boxes [44]. 
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allows one to calculate any coefficient in terms of the product of four tree-level amplitudes 
3~3|] . Up to a numerical factor, every one-loop box coefficient is given by 



B = J2AlfA^A^fA^f, (1.3) 

where the sum is over the solutions to the delta function equations and over all particles 
that can propagate in the loop. A straightforward application of this idea can be made 
for split double boxes (see fig. lb). Again the idea is to cut all eight propagators, i.e., an 
octa-cut which localizes the two loop integrations and gives the coefficient as the product 
of seven tree-level amplitudes (up to a numerical factor) 

Naively, one might expect that the coefficient of double boxes in fig. la cannot be computed 
in a similar manner. The reason is that there are only seven propagators and a hepta-cut 
does not localize the integrals over the loop momenta. 

A way to avoid the remaining integration arises in an unexpected manner. In studying 
singularities of Feynman integrals, one computes the discontinuity of an integral across a 
singularity by cutting propagators. When one cuts all propagators in a Feynman diagram 
one is computing the discontinuity across the leading singularity of the integral. However, 
at two (and higher) loops one finds a surprise when some of the external legs are massless. 
At two loops, if any of the two boxes in fig. la has at least two adjacent three-particle 
vertices (condition that is satisfied trivially for less than seven external gluons), then the 
integral has an extra propagator-like singularity beyond the naive leading singularity. The 
discontinuity across the new leading singularity is actually computed by an octa-cuti. This 
octa-cut precisely localizes the loop integrations and allows a straightforward computation 
of the coefficient as the product of six tree-level amplitudes. Up to a numerical factor, it 
is given by 

i=l 

The only two-loop amplitude in pMSYM known in the literature is the four-gluon 
amplitude J0J . One reason is that very few double scalar box integrals are known explicitly 



2 For more general double boxes, there is also an extra singularity, these are known as second- 



type singularities p5| . They cannot easily be used to produce an octa-cut but they might give a 
generalization of it. 
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HJ. In particular, to our knowledge, not all double box integrals needed for a five-gluon 
amplitude are known. Nevertheless, we present the computation of several five-gluon and 
six-gluon non-MHV scalar double box integrals as illustrations of our technique. 

This paper is organized as follows. In section 2, we review pMSYM at tree-, one-, 
two- and three-loop levels as well as the ABDK conjecture. In section 3, we show that 
the four-gluon amplitude of pMSYM can be found by using hepta-cuts. Even though the 
number of cut propagators is less than the number of integration variables, the integrand 
turns out not to depend on the loop momenta and can be pulled out of the integral. In 
section 4, we demonstrate that a certain class of double-box configurations admit an extra 
propagator-type singularity. Cutting this singularity allows us to write a universal formula 
for many double-box coefficients. In section 5, we illustrate our technique via various 
examples including non-MHV amplitudes. In section 6, we discuss applications to three- 
and higher-loop amplitudes. In particular, we show that by studying singularities, it is 
possible to realize that the basis of integrals has to contain integrals with some non-trivial 



factors in the numerator. This is in agreement with results of |44j] . 

Throughout the paper, we use the following notation and conventions along with those 
of |1| and the spinor helicity- formalism P?| , ^8| , |49[] . A external gluon labeled by i carries 



momentum Ki. Since Kf = 0, it can be written as a bispinor {Ki) a d = Xi a Xi „. Inner 
product of null vectors p a a — A a A„ and q a a = X' a X' d can be written as 2p ■ q = (A, A') [A, A'], 
where (A, A') = e a b\ a X' b and [A, A'] = eaj,A a A /& . Other useful definitions are: 

=Ki + K i+1 + ... + Kj 
K \r] = Ki + Ki+1 + ... + Ki+r _ ± 

(fl^tfrljl^EC r)[rj] 

r r 
r s r s 

mcE^xE*-)!'] = EE^ r i< r s ^ s 3\ 

r s r s 

(^(E^xE^xE^i = EEE<* r >t r s ^ s ft 

r s t r s t 

where addition of indices is always done modulo n. 



2. Review of J\f = 4 Amplitudes 

In this paper we consider amplitudes of gluons in Af = 4 super- Yang-Mills. Each 
gluon carries the following information: momentum p a „, polarization vector e a d and color 



index a. The color structure can be striped out by a color decomposition [pO| , pTl , [5^ ,p3 
Here we only consider the leading color or planar part of the amplitudes. The information 
in momentum and polarization vectors can be encoded in terms of spinors A, A and the 
helicity of the gluon h. 

2.1. Tree- Lev el Af = 4 Amplitudes 

At tree-level, the leading color approximation is exact. An amplitude is given by 

UP " e " ffi,|j yYM ^ v ' ({A<T(i),A CT (i)^ CT (i)},---,{A (T („),A (T( „ ) ,^ (T( „ ) }) 

aES n /Z n 

(2.1) 

Here we are suppressing a delta function that imposes momentum conservation. 

It is convenient to denote the set of data {Aj, Aj, hi} by i hi , where hi = ± is the helicity 
of the i th gluon. The amplitudes on the right hand side of ( |2.1|) are known as leading color 
partial amplitudes and are computed from color-ordered Feynman rules. One can study a 
given order A(l 1 , . . . ,n hn ) and the rest can be obtained by application of permutations, 
a. 

The partial amplitude A(l hl , . . . ,n hn ) can be computed using a variety of methods 
(see |S6[ for a nice review on many of the techniques developed in the 80's and 90's). More 
recently, two new techniques became available, namely, MHV diagrams 0] and the BCFW 
recursion relations [jnyn|. The latter is a set of quadratic recursion relations for on-shell 
physical partial amplitudes of gluons. For a recent review see [B4] . 



2.2. One-Loop Af = 4 Amplitudes 



Amplitudes of gluons at one-loop admit a color decomposition [ ^U| , |oT| , [5^ , |5B| , |5o| with 
single and double trace contributions. As mentioned in the introduction we will only 
concentrate on the leading color partial amplitudes 



3 It is interesting to note that since for TV = 4 SYM all particles in the loop are in the adjoint 
representation, all sub-leading color amplitudes are given as linear combinations of the planar 



ones with permutations of the gluon labels (See section 7 of [42] for a proof. 
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One-loop amplitudes of gluons in supersymmetric theories are four-dimensional cut- 
constructible p2| , |56|1 . This means that they can be completely determined by their finite 
branch cuts and discontinuities. N = 4 amplitudes are even more special. Reduction 
techniques ]43[ can be used to express these amplitudes in terms of scalar box integrals 



||42|| . These are one- loop box Feynman integrals in a scalar field theory where a massless 
scalar runs in the loop, 

I = I dH {i 2 + ie)({£ - ki) 2 + ie)((£ - k[ - k 2 ) 2 + ie)((£ + k 4 ) 2 + ie) (2 ' 2) 

where ki, k 2 , k%, k 4 are the external momenta at each vertex. They are not independent 
since by momentum conservation /C3 = — (£4 + k% + k 2 ). Note that the integral ( |2.2[ ) is 
singular when at least one ki is a null vector. Therefore, we should specify a regularization 
procedure, like dimensional regularization. However, we will be considering cuts that 
are finite and do not depend on the regularization procedure. Since A(l, . . . , n) is color- 
ordered, each k can only be the sum of consecutive momenta of external gluons. Moreover, 
since we only consider the planar contributions, we can define a given contribution by 
specifying k, I such that k\ = Ki + . . . + Kj-i, k 2 = Kj + . . . + Kk-i and ks = 
Kj~ + . . . + The reduction procedure then gives for the amplitude an expansion of 



the form 42 



A(l,...,n)- 2J B ijkl Z(K i +...+K^ 1 ,K :i +...+K k _ 1 ,K k +...+K l _ 1 ), (2-3) 

Ki<j<t<m<n 

where the coefficients Bijj-m are rational functions of the spinor products. Since all scalar 
box integrals are known explicitly, the problem of computing A(l, . . . , n) is reduced to that 
of computing the coefficients B^i. 

A general formula for the coefficients was found in in terms of products of 



tree level amplitudes. Let us review the derivation of the formula because the idea is useful 
in the analysis at higher loops. If we think of the scalar box integrals as an independent 
basis^ of some vector space we can interpret A(l, . . . , n) as a general vector. All we need 
to do is to find a way to project A(l, . . . , n) along the space spanned by a given scalar box 
integral I. From the definition of X in (|2.2|) we see that each integral is uniquely determined 
once its four propagators are given. It is natural to think that the way to determine the 
coefficient B is by looking at the region of integration where all four propagators become 



The notion of independence is the equivalent of cut constructibility of the amplitude. 
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singular. In fact, the integral obtained by cutting, i.e., by dropping the principal part of 
all four propagators computes the discontinuity of the given scalar box integral across a 
singularity which is unique to it. 

The set of four equations that gives I is the following 



0, 



o, 



(£-h-k 2 f 



0, 



+ kA 2 



0. 



(2.4) 



A little exercise shows that these equations do not have a solution if I is a real vector 
in Minkowski space for general external momenta. The way out of this problem is to 

complexify all momenta and make a Wick rotation to ( h+) signature. In the new 

signature the delta functions are still well defined and there are always solutions to (|2.4|) . 




m+1 



Fig. 2: A quadruple cut diagram. Momenta in the cut propagators flows clockwise 
and external momenta are taken outgoing. The tree-level amplitude Ai ree , for 
example, has external momenta {i + 1, j, £2, ii}- 



One can also look at the same regime on the left hand side of ( |2.3| ) by considering only 
Feynman diagrams that posses the four propagators entering in ( ^.4j ) . Summing over them 
one finds the following equation 

J dfi J2 njA^fA^fA^fA^f = B ijkm J dfi (2.5) 

where sum over J represents a sum over all possible particles in the M = 4 multiplet. The 
measure dfi is the same one both sides of the integrals, 

dv = d 4 £ 5 {+) {f) 5 {+ \{l- h) 2 ) S^di-h-h) 2 ) 5 i+) ((£ + k 4 ) 2 ), (2.6) 



and the tree-level amplitudes are denned as follows (see fig. 2) 

A tree j + 1, i + 2, . . . , J - 1, J, € 2 ), Agf = A(-£ 2 J + 1, J + 2, . . . , k - 1, k,£ 3 ), 

A%f =A(-£ 3 , k + l,k + 2,...,m-l,m,£ 4 ), A\% e = A(-£ 4 , m + 1, m + 2, . . . , i - 1, i, £ x ) 

(2.7) 

where 

*i=t, h=£-kx, £ 3 = £ ~ h - k 2 , h = £ + k 4 , h = K l+1 + . .. + Kj, 

(2.8) 

k 2 = K j+1 + . . . + K k , k 3 = K k+1 + . . . + K m , k 4 = K m+1 + . . . + K t . 

The integral J d/u, is just given by a Jacobian 1/vA. This Jacobian cancels on both sides 
since the integral is localized by the delta functions and the coefficient is given by 

B ijkl = ±J2 n J A t?) eA lW A %T A T4 e - ( 2 - 9 ) 



S,J 



Here S is the set of solutions to the conditions imposed by the delta functions, and |<S| 
is the number of solutions. The sum also involves a sum over all possible particles that 
can propagate in the loop. For further details and many examples we refer to ]33|]. Even 
though the Jacobian did not play an important role for the quadruple cut technique at 
one-loop, it is crucial for the two- loop analysis we carry out in section 4.1. For this reason 
let us write it down for future reference 

A = s 2 t 2 - 2st{k{kl + k\k\) + (kfkj - k\k\f (2.10) 

with s = (k\ + k 2 ) 2 and t = (k 2 + k 3 ) 2 . 

2.3. Two-Loop Af = 4 Amplitudes 



At two loops, only the four-gluon amplitude has been computed |[44[| . The M = 4 
calculation was the first full two-loop amplitude of gluons ever computed. The answer 
is given as a linear combination of double box scalar integrals with coefficients that are 
rational function of the spinor variables. A double box scalar integral is the analog of the 
one-loop scalar box integral introduced above, more explicitly, 

f _d*p 1^ 

(A*,..-, 6)-J (27r )4 (p2 + ie)((p _ kl) 2 + ze)((;p _ ^ _ k2) 2 + lt) X 

A i^ 

(2tt) 4 ((p + q + k 6 ) 2 + ie)(q 2 + ie)((q - k 5 ) 2 + ie)((q -k±- k 5 ) 2 + it) ' 

(2.11) 
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This integral is UV finite but it might have IR divergences when some fc's are null vectors. 
Again, as in the one-loop case, one has to choose a regularization procedure but we do 
not do so because we only discuss finite cuts. The planar contribution to the four-gluon 



amplitude is 



A 2 - loop (K 1: K 2 , K 3 , K 4 ) = A\ ree s t (s X{K U K 2 , 0, K 3 , K 4 , 0) + t I(K 4 , K u 0, K 2 , K 3 , 0)) 

(2.12) 

where s = (K l +K 2 ) 2 and t = (K 2 +K 3 ) 2 . This was computed by using the unitarity-based 
method [fil^|5l||5"T|,|5*£|,|5T)|| . It is very important to mention that the double box scalar integral 



( [2. 11|) is not known in general but explicitly formulas exists in dimensional regularization 



when k 3 = k§ = and at least three of the other kiS are null vectors [46 



2.4. ABDK Conjecture 

As mentioned in the introduction, one of the motivations of this work is to prepare 
the ground for a more extensive test of the ABDK conjecture. This conjecture asserts that 
the planar limit of L-loop amplitudes in M = 4 SYM is determined iteratively, i.e., as a 
function of /-loop amplitudes with I < L. 



Let us make this more precise. Here we follow |37| and [[ID] where the original proposal 
was made. Consider the function 

M ^ *= A ^£:$ ^ 

then the ABDK conjecture states that 

M( i '=P i (MW ) ...,M( i - 1 ') (2.14) 

where Pl(xi, ■ ■ ■ , ^l-i) is a certain polynomial of degree L and independent of the helicity 
configuration. The explicit form of ( |2.14j ) at two loops was given in in terms of the 
function /(e) = (-0(1 — e) — ip(l))/e, where the digamma function is defined by ij){x) = 
r'(x)/r(x), as follows 

M?\e) = \ (Mi 1 \e))\ /(e)M«(2e) - |c*. (2.15) 

This conjecture was explicitly checked for four-gluon amplitudes. Very recently, the form 
of the polynomial in ( |2.14| ) was obtained for the three-loop four-gluon amplitude in |4l[ . 
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One of the impressive predictions of the conjecture is a relation between the finite re- 
mainders which are defined at e = 0. At two loops, one introduces the universal singular 
function C n (e)^ [ pDIJS 7 ^ ] which contains the infrared singularities and does not depend on 
the helicity configuration since it is normalized by the tree-level amplitude. Defining the 
finite remainder as 

Fi 2 \e)=Mi 2 \e)-C^(e), (2.16) 



one can write a finite (as e — > ) analog of (|2.15|) as follows [[37 

'lira 



4 2) (o) = ^(Fi 1) (o)) 2 -C 2 4 1) (o)-^ 



+ 5 C 



(2.17) 



Recall that at one- loop F^\0) can have at most dilogarithms, while F^\o) can have 
higher polylogarithms. This means that very non-trivial cancelations must happen. These 
cancelations were found to occur for n = 4 between terms coming from the two integrals 
in (|2.12|) and involved many polylogarithmic identities [37j. In the recent paper |[41|| , an 
impressive formula for the all loop finite remainder of MHV amplitudes was also presented. 
The formula is given in a kind of generating function structure 



l + ^a L F,( L) (0)=exp 



L=l 



l>y K F<p(0)+C 



(2.18) 



where a is basically the 't Hooft coupling, 'Jk is the universal soft anomalous dimension 
and C is a function that admits a power series representation in a. 'Jk and C are known 
up to the order needed to obtain the three-loop termi. 



3. Four-Gluon Two-Loop Amplitudes and Hepta-Cuts 

In this section, we consider hepta-cuts of the two-loop four-gluon leading partial am- 
plitude. This section can be viewed as a warm-up section where we introduce relevant 
notation and do some sample calculations which will be used in the rest of the paper. It 
is enough to consider Al~ loop (l~, 2~, 3 + , 4+) as all other A 2 4 ~ loop with different helicity 



5 It is important to mention that there is no canonical definition of the finite remainder F n . 
In fact, the definition of finite remainder used in ( [2.17| ) differs from that used in ( [2.18| ). For more 



details see [41]. We thank Z. Bern and L. Dixon for useful discussions on this point. 
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assignments can be obtained from this one by Ward identities. The leading partial ampli- 



tude was first computed in pH}] . As reviewed in section 2.3, the amplitude can be presented 
as a linear combination of two scalar double-box integrals (see fig. la) 



J 



d 4 p d 4 q 



(2tt) 4 (2tt) 4 p 2 (p - Krfip -K x - K 2 ) 2 (p + q) 2 q 2 (q - K 4 ) 2 (q - K 3 - K±) 2 



(3.1) 



where Ki are the four external gluon momenta, with rational coefficients. All external 
momenta are assumed to be outgoing. The integral ( |3.1| ) has seven propagators, hence 
it is natural to consider hepta-cuts. It turns out that the coefficients can easily be found 

from hepta-cuts when the loop momenta are analytically continued to signature ( h+) 

or complexified. In the present case, there are two independent coefficients as well as two 
independent hepta-cuts. We refer to them as the s-channel cut and the t-channel cut. The 
corresponding coefficients will be denoted as c s and q. Let us start with the cut in the 
s-channel. In this case, there are six different helicity configurations. For all of them, only 
gluons can propagate in both loops. A sample helicity configuration is shown in fig. 3. In 
this paper, for simplicity, we depict tree level amplitudes as points. Since all propagator 
are cut, there is no need to indicate a cut by a dash line and we choose to omit them i. 



K 4 .+ 



■K, 



_ q 



+ - 



+ 

<Jiv 



Al 



+ 

/Pi 



+ 



q 2 



+ 



p 2 



+ 



K 



+/ 



^3 ,v 2 
Fig. 3: A sample hepta-cut in the s-channel. Tree level amplitudes are depicted 
as points. All propagators are cut and therefore we omit the dashed lines used in 
fig. 2. 



Note that conventions in fig. 3 are different from those used in fig. 2 where all tree level 
amplitudes are denoted by blobs and cuts are indicated by dashed lines. 
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The rational coefficient c s is then given by 



7 f A,, I Atree Atree Atree Atree Atree Atree\ 

2^1=1 J "AH^m ^(2) ^(3) A U) ^(5) ^(6) )i 



^-x„ . ; (1) ^(2) ^(3) ^(4) ( 6 ) - /o 9 ^ 

where by A f ^ e we denote tree amplitudes at each of the six vertices, the integration 
measure dfi is given by 

dn =^-0^ 6 (P 2 W(P - KifWiP ~ K i- K2) 2 )S((p + qf) 
5(q 2 )5((q - K 4 ) 2 )5((q - K 3 - K 4 f) 

and in the numerator we sum over the product of the six tree level amplitudes corresponding 
to a given helicity configuration. The factor i r comes from the seven propagators. It seems 
that since the number of delta functions is less than the number of integration variables, 
the integral does not localize and one integration has to be performed. However, it turns 
out that the integrand can be simplified in such a way that the dependence on the loop 
momenta drops out and we are left with the integral of the measure which cancels out 
according to eq. (|3.3|). 

In the discussion of one-loop amplitudes in section 2.2, we mentioned that the mo- 
mentum I has to be complexified in order to find solutions to the four equations from the 
cut propagators. Making I complex also has as a byproduct the fact that three-particle 
vertices on-shell do not have to vanish. Tree-level three-gluon amplitudes with helicities 



[ h) or (+ H — ) are given respectively by [|61| , |62 



A t r £ (p-,q-,r+)=t {P f Ar e (p + ,q + ,r-) = -i T ^— y (3.4) 

\q r) {r p) [q r\ [r p\ 

In Minkowski space, A p and A p are related to each other as A p = ±A P . This means that if 
p ■ q = 0, which follows from momentum conservation at the vertex, then both (A p X q ) = 
and [A p X q ] = 0. This implies that both amplitudes in (|3.4|) vanish. If we complexify the 
momenta, then the equation p ■ q = has two independent solutions. We have that either 
(A p X q ) = or [Ap X q ] = 0. That is either A p and X q are proportional or A p and X q are 
proportional. Also note that momentum conservation implies that p-q = p- r = q- r = 
0. This means that either three A's are proportional or all three A's are proportional. 
Therefore, for every (+ H — ) tree level amplitude we choose all A's to be proportional. 
Similarly, for every ( h) tree level amplitude we choose all three A's to be proportional. 
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Explicit calculations, considered for one of the helicity configurations in some detail 
below, show that every helicity configuration gives the same contribution equal to 

-Al ree s 2 t J dp, (3.5) 

where A l l ee is the tree-level four-gluon amplitude 

A* ree (l-,2",3 + ,4+) = i- -, (3.6) 

4 v ' ' ' ; (2 3) (3 4) (4 1)' v ; 

and 

s = (K x + K 2 ) 2 , t = (K 2 + K 3 ) 2 . (3.7) 

Note that the integral in ( |3.5| ) cancels against the denominator in ( |3.2| ) . The coefficient 
6 in the numerator will also cancel. The reason is the following. In the denominator in ( |3.2| ) , 
we have to sum over all different solutions to the delta-function conditions. It is easy to 
realize that in this particular case the number of different solutions equals the number of 
helicity configurations. Thus, each term in the numerator in ( |3.2|) picks one of the six 
solutions whereas in the denominator we sum over all the six solutions. As a result, we 
obtain 

Cs = -A\ ree s 2 t. (3.8) 



This coincides with the corresponding coefficient found in 

Let us consider the helicity configuration shown in fig. 3 in some detail. The analysis 
of the remaining five configurations is completely analogous. Consider the integrand as 
the product of six tree-level amplitudes 



[PUP? \ ( (Pi,2) 3 

h i][i,piL 

fe,/] 3 \ ( (pJ) S \ ( [9i,4] 3 \ / (q u q 2 ) 3 



•I a tree a tree a tree a tree a tree Atree\ 

A m A (3) A W A {5) A {6) ) 1 _-» v M[ljpi]j V(2 , P2 )(p 2 , Pl ) 



(3.9) 



[l,P2}\p2,q2}J \(l,q)(q,p)J \[^q}[Q:Qi}J \ (92, 3) (3, ql) 

Next, simplify this expression by using momentum conservation. For example, the product 
of \pi,p] and (pi, 2) can be simplified as follows 

[pi,p](pi,2> = -{2\ Pl \p} = (21^] = -(1 2)[1 p]. (3.10) 

Then the product of the first four factors in (|3J]) becomes 

(12} 2 [q,q 2 ] 2 . (3.11) 
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After using momentum conservation along the lines of eq. (|3.10 ), one finds 



•/ a tree Atree Atree Atree Atree Atree\ ■ 2 (! 2 ) 2 [ 3 4 ] Atree 2+ /o in^ 

_ H A (1) ^(2) A (3) A (4) A (5) A (6) )l — lS — ~ A 4 s l - (o.l-J 



Note that this expression does not depend on the loop momenta and, thus, can be pulled 
out of the integration. 

Now we consider the hepta-cut in the t-channel. Here we have ten helicity configu- 
rations. Note that in this case the number of helicity configurations does not equal the 
number of solutions of the delta-function equations. By a solution we mean a choice 
whether all A's are proportional or A's are proportional at each of six three gluon vertices. 
However, among the ten configurations, there are different configurations for which the 
choices of whether A's or A's are proportional are exactly the same. A solution then means 
summing up over such configurations. In this case, there are two helicity configurations 
corresponding to actual solutions and the remaining eight ones break up in pairs. The sum 
of the two helicity configurations in each pair corresponds to a solution of the delta-function 
conditions. Overall, we have six improved helicity configurations, each corresponding to an 
independent solution to the delta-function conditions. All paired up configurations involve 
fermions and scalars running in one of the loops and summation over the two configurations 
in a given pair provides a significant simplification. The coefficient Ct is given by 

„■ V^IO r j ( Atree Atree Atree Atree Atree Atree\ 
_ 1 2^1=1 J a H\ A (i) A {2) ^(3) A (4) ^(5) ^(6) I 1 f0 10 v 

Ct ~ J> ' (3 - 13) 

where the sum is over all ten configurations, or over the six improved configurations, each 
corresponding to an actual solution of the delta-function equations. As before, all six 
improved configurations give the same contribution 

-A^st 2 . (3.14) 

Since each improved configuration corresponds to a solution to the delta-function equa- 
tions, the factor 6 cancels out. As a result we obtain 

^ = -Al ree st 2 . (3.15) 

This coincides with the corresponding coefficient from [ 13 1 . As an example, let us consider 
the two helicity configurations shown in fig. 4. 
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Fig. 4: Examples of hepta-cuts in the i-channel that correspond to the same 
solution of the delta function constraints. 



Note that for both double boxes the choices whether all A's or all A's are proportional 
are exactly the same at every vertex. Therefore, it is the sum of these two diagrams that 
corresponds to one of the six improved helicity configurations. Both double boxes in fig. 4 
involve gluons, fermions and scalars running in one of the loops and only gluons running 
in the remaining loop. The necessary tree level amplitudes are given by 

A^(p-, g- r+) = i {pq)3 ( {q V) ) a , A tree (p+, q+, r~) = -1 J^l!_ ( M 
" ' (q r){r p) \(q p) J ' ' [q r] [r p] \ [q p] 

(3.16) 

where a = for gluons, a = 1 for fermions and a = 2 for scalars. After summation over 
the two configurations, we obtain 



-i , (3.17) 

7 



where 



a = (qi 1) (92 P2) b q] [4 ?! , P = (1 g 2 ) (P2 / P g 4] 

(3.18) 

7 = <?2 gi)(i ga>[5 i][gi g](i ?2) (P2 p][g 4](g l)(g 2 p 2 )[4 gi][p gj. 

By using momentum conservation along the lines of eq. ( |3.10|) , we can simplify ( |3.18| ) 
to obtain —A tree st 2 . The integral of the measure factors out and cancels against the 
denominator in (|3.13| ) to give ( |3.15| ). 

Thus, we find that the coefficients of the four-gluon amplitude double boxes can be 
calculated by studying hepta-cuts. Of course, this is not enough to claim that this is the 
full answer. One still has either to prove that the answer has all the correct discontinuities 



across all branch cuts, which was done in [44], or to prove that the basis of integrals is 



15 



given entirely by double boxes. Unfortunately, the basis of integrals is not known at two 
loops. 

For four gluons, even though the number of integration variables is greater that the 
number of the delta- functions in a hepta-cut, no integration has to be performed. We find 
that this is not the case if the number of external gluons is greater than four. We will 
see in the next section that already in the case of five-gluon amplitude, the product of 
the corresponding tree-level amplitudes does depend on the loop momenta and cannot be 
pulled out of the integral. 

4. Octa-Cuts of Two-Loop Amplitudes 

In the introduction we distinguished between two different kinds of double box scalar 
integrals. In the first class, the two boxes share a propagator while in the second class 
they only share a vertex, see fig. la and fig. lb respectively. In this section, we show how 
one can use octa-cuts to compute the coefficient of a certain subset of the first class and 
the coefficients of all integrals of the second class, which we called split double boxes. 

4-1. Double- Box Scalar Integrals 

Let us start with the double boxes that have seven propagators. We will show that 
when at least one of the two boxes has two adjacent three particle vertices then there 
is an extra propagator-like singularity that can be cut. This produces one more delta- 
function which together with the hepta-cut of the previous section completely localizes 
the cut integral. Even though we concentrate only on the planar configurations, exactly 
the same logic can be applied for non-planar configurations as well. Consider an arbitrary 
double-box configuration shown in fig. 5. The corresponding hepta-cut integral is 

1 = J ( ^4 ^ (p 2 mp-ki?mp-^^ 

(4.1) 
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Fig. 5: An arbitrary double-box configuration. 



Let us perform, for example, the p-integration. The integral over p, 

X p = J d* P 5(p 2 )5((p - h) 2 )5((p -h- k 2 ) 2 )5((p + q + k G 



is localized and the answer is 



(4.2) 



T 



where 



P = 
Ai 



/ 1-2(A 1 +A 2 ) + (A 1 -A 2 ) 2 
k 2 (k 3 + k 4 + k 5 - q) 2 



{k 1 + k 2 Y(k l + k & + q y 
The crucial observation is that when 



A, 



k 2 (k 6 + qf 



(h + hnh + h + q)*' 



(4.3) 



(4.4) 



(4.5) 



X acquires an extra propagator-type singularity, i.e. l/(&i + ke + q) 2 ■ We can formally 
cut the new propagator by replacing it with a delta-function creating an eighth cut. In 
other words, after performing the p-integration we end up with following integral over q 
(we omit the overall g-independent factor) 



(4.6) 



l q = [ d 4 qS(q 2 )S((q - k 4 ) 2 )5((q - k 3 - k 4 ) 2 )— — i 
J («i + h 

This integral looks like a triple cut of the following effective box 



6 + q) 2 
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Fig. 6: Effective box that arises after a quadrupie cut is used to iocaiize the p 
integral. The momentum flowing along the uncut line is q + k\ + k%. 



Note that the momentum flowing along the uncut line is exactly q + k\ + k$. From this 
viewpoint it is natural to cut the remaining propagator. Note that this procedure localizes 
the g-integral. Then it is straightforward to write down the coefficients of such double-box 
integrals. They are given by 

r — % _ \ (n t n r Atree Atree Atree Atree Atree Atree\ ( A 7\ 

a ~ |£| \ n Ji n J2 A {l) A {2) ^(3) ^(4) ^(5) ^(6) )hi K^-l) 

h,Ji , ^2i<S 

where the sum over h is the sum over all helicity configurations, the sums over J\ and Ji 
are the sums over all particles that can propagate in both loops, S is the set of all solutions 
for the internal lines of the following system of equations 

p 2 =0, (p-kx) 2 = 0, (p - h - k 2 ) 2 = 0, (p + q + k G ) 2 = 0, 

(4.8) 

q 2 =0, (q-k 5 ) 2 = 0, (q - k 4 - k 5 ) 2 = 0, (£q + k 6 + q) 2 = 0, 
and \S\ is the number of solutions. This expression is analogous to the formula for one- loop 



coefficients of box integrals [|34[| . It is important to remember that this discussion is valid 
if 

p = l, Ai = A 2 = 0. (4.9) 

Otherwise, the singularity l/(fci + k§ + q) 2 will be replaced by a more complicated one 
which is not propagator-like, as it can easily be seen from eq. (|4.4|) . The conditions given 
in ( f4.9| ) are satisfied if a given box has two adjacent three-particle vertices. It easy to check 
that this is always the case if the number of gluons is less than seven. This means that 
every double-box coefficient of any gluon amplitude with less than seven external lines 



is given by eq. (|4.7|) . The first double-box configuration where eq. Q4.9p is not satisfied 
appears when the number of external gluons is seven and is shown in fig. 7. 
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Fig. 7: The simplest double-box configuration for which the conditions in ( f4.9| ) 
are not satisfied. 



However, even if the number of external gluons is greater than six, there are double-box 
configurations for which eq. Q4.9] ) is satisfied. In such cases the eighth cut still exists and 
eq. ( |4.7|) is still correct. 



In fact, eq. (|4.7| ) requires some additional explanations. Note that existence of the 
effective box in fig. 6 implies that either the momentum I or the momentum pi in fig. 
5 vanishes. In Minkowski space, this would mean that some tree level amplitudes in 



eq. (|4.7| ) vanish. Moreover, in Minkowski space, the system of equations ( |4.8| ) does not 



have solutions, which means that we cannot see the singularities under consideration. 
Therefore, it is not surprising that eq. ( ft.7| ), at least naively, is meaningless in Minkowski 
space. In order to see the new kind of singularities, we have to analytically continue 

all momenta to signature ( 1"+)- But in signature ( the statement that a 

tree amplitude vanishes when one of the incoming or outgoing momentum vanishes is not 
correct. Each tree amplitude is constructed by using spinors. When one of the incoming or 

outgoing ( h+) momenta vanishes, it is impossible to determine its spinors components 

even up to rescaling. This leaves the amplitude undetermined. For example, assume that 
we have a helicity configuration containing a three-gluon amplitude A(p~ 1 p^[ 1 k^~). It is 
given by 

(4.10) 



(pfci)(fciPi) 



If pi vanishes, the spinor A Pl cannot be uniquely determined. In fact, A Pl is not uniquely 
defined even for non-zero pi as it is defined up to rescaling. However, when pi = the 
freedom in not being able to determine A Pl becomes much larger. One can always say 
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that pi = implies that X Pl = and A Pl is arbitrary. This means that A(p~ 1 pi,kf) 
becomes arbitrary. Therefore, the numerator in eq. ( [4.7| ) is a discontinuous function of 
momenta and we have to give a prescription on how to define it as / or pi goes to zero. 
The natural way to define it is as follows. Consider first the loop with momentum p. Let 
^4*i) e , ^.(2) e ? ^(3) 6 an d ^-(4) 6 be the four tree amplitudes which depend on p. Assuming 
that they are all non-zero, we can solve the first four p-dependent equations in ( |4.8|) to 
determine p as a function of the external momenta and q and then evaluate the product 
A l ^ e A%^ e A^ e A^ e on these solutions. We claim that this product can be simplified in 
such a way that it is a well-defined function when the constraint {k\ + k$ + q) 2 = is 
imposed. Below, we will present a few examples that show that this is indeed the case. 
Having found the product A l ^ e A l ^ e A l ^ e A 1 ^ as a function of the external momenta and 
q, we then multiply it by the remaining two tree amplitudes A^ e and A^ e and evaluate 
the product on the solution of the remaining four equations in (f4.8|). We propose this as 
a method for calculating double-box coefficients provided conditions ( }L9| ) are fulfilled. 

A Subtlety 

There is one important subtlety we have to discuss before presenting examples. Con- 
sider a helicity configuration with two adjacent three-particle vertices, one of which depends 
only on internal momenta and the other one depends on both internal and external mo- 
menta, with both vertices having the same helicity configuration. For example, consider 
the configuration in fig. 8. 
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This configuration is non-zero only if \ q ~ Ai. Therefore, the integral over p 




tree 



(4.11) 



must be proportional to 8{{k\ + q) 2 )- In other words, the integral lacks the extra 
propagator-like singularity and therefore it does not contribute to the octa-cut. 

4-2. Split Double-Box Scalar Integrals 

When the number of gluons is greater than five, a new kind of double box integrals can 
appear. These were introduced in the introduction in fig. lb. For the reader's convenience 
we depict them again in fig. 9. This double box scalar integrals are such that the two 
boxes only share a vertex and not a propagator. This is why we will call them split double 
boxes. 



Fig. 9: Generic split double box configuration. 

The coefficients of the split double boxes are easy to compute. Since they have eight 
propagators and the two loop integrations are completely independent, it is straightforward 
to consider two quadruple cuts or equivalently an honest octa-cut. This produces eight 
delta-functions that localize both loop integrals. Let us see this in more detail. The 
quadruple cut in the g-loop fixes q to be a solution to the following equations, 




q 2 = 0, (q- k 6 ) 2 = 0, (q-k 5 - k 6 ) 2 = 0, (q + k 7 ) 2 = 0, 



(4.12) 



while a quadruple cut in the p-loop fixes p to be a solution of 



p 2 = 0, (p-k 2 ) 2 = 0, (p-k!-^) 2 = 0, (p + k 3 ) 2 = 0. 



(4.13) 
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Each set of equations has two solutions. The coefficient of a split double-box scalar integral 
is then given by 

r—— (n 7 n t Atree Atree Atree Atree Atree Atree Atree\ (a -ia\ 

h,Ji , J2 ,S 

where we have used that the number of solutions is 4. 
5. Examples 

In this section we consider several examples of coefficients calculated by using octa- 
cuts. All of them are coefficients of scalar double boxes with seven propagators. We 
consider four-, five- and six-gluon amplitudes. For six-gluons we study a non-MHV ampli- 
tude with adjacent negative helicity gluons. 

5.1. Four- Gluon Amplitude A 2 ~ loop (—,—,+, +) Revisited 

As a first example, let us reconsider the octa-cut of the four-gluon amplitude from 
section 3 in the s-channel. The octa cut in the t-channel is analogous. The additional 
propagator that we cut is j-j^^p. . See fig. 3 for notation. Taking into account the 
subtlety in the previous section, there are only four helicity configurations that contribute. 
All of them give the same answer — A t £ ee s 2 t. On the other hand, the number of solutions 
to eqs. (|4.8| ) can be shown to be four. This gives c s = —A t l [ ee s 2 t as in eq. (|3.8|) . Note that 
the product of the four tree level amplitudes depending on the internal momentum p is 
given by (1 2) 2 [q q^] 2 (see eq. ( |3.11| )). This expression does not have any ambiguity in the 
presence of the eighth delta- function d((q + Ki) 2 ). 

5.2. Five-Gluon Amplitude A 2 ~ loop (-, -, +, +, +) 

As a next example, let us calculate the coefficient of the following five-gluon double- 
box configuration. 

In this case, there are two helicity configurations to consider. Both of them can be shown 
to give the same contribution. We will consider the helicity configuration shown in fig. 
11. 
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Fig. 10: A double box integral of the five-gluon amplitude A(l ,2 , 3 + , 4 + , 5 
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Fig. 11: One of the two possible helicity configurations contributing to the coef- 
ficient of the integral of fig. f 0. 



Note that only gluons can propagate in both loops. The product of the six tree level 
amplitudes is as follows 



(1 2) 2 [g q 2 ]' 



[3 qi] 



(q QiY 



(5.1) 



[qi q2][q 2 3] ( Ql 4) (4 5) (5 q) 
where the first two factors come from the four vertices on the right. The computation 



leading to the first two factors was done in the previous section in eq. (|3.11| ). By using 
momentum conservation along the lines of eq. ( |3.10|) , eq. (|5.1|) can be reduced to 

(1 2) V 



(3 4) (4 5) 



[5 3] + [4 3] 



(g 4) 
(q 5) 



(5.2) 



Now we impose the delta-function S((Ki + q) 2 ). It has two possible solutions, X q ~ Ai or 
A„ ~ Ai. It is not hard to show that if we choose A„ ~ Ai then the expression in eq. 
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vanishes. Therefore, the only relevant solution is A„ ~ Ai. Then the octa-cut becomes 



(3 4) (4 5) \} b 3J + [4 3J (15) 



d\i. 



(5.3) 



Taking into account that the system of equations (|4.8| ) has four solutions, we find that the 
corresponding coefficient is 

2% (1 2) 2 s 2 



,(5) 



4 (3 4) (4 5) 



[5 3] + [4 3] 



(1 4) 
(15) 



1 



--At ree s 2 t 



(5.4) 



where s = (K\ + K2) 2 and t = (K 2 + K3) 2 and A^~ ee is the tree-level five-gluon amplitude 

(1 2) 3 



Af ee = i- 



(5.5) 



(2 3)(3 4>(4 5)(5 1) ' 

Let us consider one more example. Let us calculate the coefficient of the five-gluon 
double-box configuration shown in fig. 12. 



+ 



-K, 



K 3 + - 



Fig. 12: Second example of a double-box configuration in the five-gluon amplitude 

A(l-,2-,3+,4+,5 + ). 



There are two helicity configurations to consider. However, one of them gives the zero 
answer. The non-zero contribution comes from the helicity configuration shown in fig. 
13. 

Note that only gluons can propagate in both loops. The number of solutions to eqs. (|4.8| ) 
can be shown to be two. Then the corresponding coefficient is given by 

,(5) _ _ i (qi Q2) 3 [qi 4] 3 (p if [q 2 I} 3 

2 (q 2 3) (3 Ql ) [4 q][q Ql ] (I q)(q 5) (5 p) [l p 2 }[p 2 q 2 ] 
{px 2) 3 [ Pl pf 



c 2 



(2 p 2 )(P2 Pi)[pl][2 Pi]' 
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Fig. 13: The only non-vanishing helicity configuration contributing to the coeffi- 
cient of the double box integral of fig. 12. 



Using momentum conservation and the fact that X q ~ A4, eq ( |5.6| ) can be simplified to 
give 

4 5) = --Af^stu, (5.7) 

where u = (K 3 + K4) 2 . All other double-box coefficients of the five-gluon amplitude can 
be found by analogous calculations. 

5.3. Six-Gluon Amplitude A 2 - loop (-, -, -, +, +, +) 

As our next example, let us calculate the coefficient of the six-gluon next-to-MHV 
double-box configuration shown in fig. 14. The additional singularity that we cut is again 
/_^ K \2 ■ There are two helicity configurations to consider, both giving the same answer. 
In both of them only gluons can propagate in both loops. Let us describe the calculation 
of the one depicted in fig. 15. 




Fig. 14: A double-box integral of the six-gluon non-MHV amplitude A(l ,2 ,3 ,4 + ,5 + ,6 + ). 
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Fig. 15: One of the two helicity configurations contributing to the coefficient of 
the scalar double box of fig. 14. 



The product of the four p-dependent tree amplitudes gives (see eq. (|3.11|) ) 

(12) 2 [^ 2 ] 2 . (5.8) 



Then the numerator of eq. (4/7) becomes 



n-( Atree Atree Atree Atree Atree Atree\ 

- A {2) ^(3) ^(4) ^(5) ^(6) > — 

(1 2)*\a a I 2 ^ 3)3 [5 ^ (5 - 9) 

{ } [qQ21 (3 4}(4g 1 )(g 1 g 2 )[6g][ ? g 1 ][ gi 5]' 

where the factor of two comes from two helicity configurations. Using momentum conser- 
vation similar to eq. ( |3.11| ), eq. ( |5.9| ) can be simplified as follows 

(12) 2 [5 6] 3 (3|5 + 6|q] 2 (q 2 3) 
(3 4>[6g](4|5 + 6|g](ga|3 + 4|5]" 1 ' ' 

Now we impose the last condition (q + Ki) 2 = 0. This equation has two solutions. We 
can either have X q ~ Ai or X q ~ Ai. Both solutions give non-zero contributions. The first 
solution yields 

. (12) 2 [5 6] 3 (3|5 + 6|1] 2 (2 3) 

(3 4) [6 1](4|5 + 6|1](2|3 + 4|5] 1 ' J 

while the second solution yields 

, (1 2) 2 [5 6] 3 (5 6) [4 2] 2 (3|(1 + 2) ■ (5 + 6) ■ (3 + 4)|2] 
[23](5|6 + l|2][5|(3 + 4).(l + 2).(5 + 6).(3 + 4)|2]" l ' J 
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Taking into account that the system (|4.8| ) in this case has four solutions, the double-box 
coefficient becomes 



,(6) 



(1 2) 2 [5 6] 3 (3|5 + 6|1] 2 (2 3) 
(3 4) [6 1](4|5 + 6|1](2|3 + 4|5] + 

(1 2) 2 [5 6] 3 (5 6) [4 2] 2 (3|(1 + 2) ■ (5 + 6) • (3 + 4)|2] \ 
[2 3] (5|6 + 1|2] [5|(3 + 4) • (1 + 2) • (5 + 6) • (3 + 4)|2] J ' 



(5.13) 



Let us consider one more example. Let us calculate the coefficient of the six-gluon 
double-box configuration shown in fig. 16. 
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Fig. 16: Second example of a six-gluon double-box integral of the six-gluon non- 
MHV amplitude A(l~, 2~ , 3~, 4+, 5 + , 6+). 
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Fig. 17: The only non- vanishing helicity configuration contributing to the coeffi- 
cient of the scalar double-box integral of fig. 16. 



In this case, there is only one helicity configuration contributing to the octa-cut. It is 
shown in fig. 17. Only gluons can propagate in both loops. The system of equations ( |Q| ) 
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has two solutions. Therefore, the corresponding double-box coefficient is given by 



.(6) 



[pi pY 



(Pi 2} 5 



(p iy 



2 [p! 1][1 Pl ] (2 P2 )( P2 Pl ) (/ p) [/ P2 ][ P2 3] [3 g 2 ] 
fei 6] 3 (qi q 2 ) 3 



x 



(5-14) 



[6 q][q qi] (q 2 4) (4 5) (5 Ql ) ' 
By using the first seven equations in (|4.8| ), we can simplify ( |5.14| ) as follows 



J6) 



u 3 s(l\q\6] 



where 



2[1 2] [2 3] (4 5) (5 6><4|g - Kf J |3] 



[K 4 + K 5 + K 6 ) 2 , s = (K 1 + K 2 ) 2 . 



(5.15) 



(5.16) 



Now we consider the last equation (q + Ki) 2 = 0. From fig. 17, it follows that X q has 
to be proportional to Xq. Therefore, X q has to be proportional to X\. Using momentum 
conservation, one can find that 

u 



(6|4 + 5|1 

Substituting this into eq. ( |5.16|) , we obtain 



■AeAi. 



(5.17) 



,(6) 



u 4 st 



(5.18) 
(5.19) 



2 2 [1 2] [2 3] (4 5)(5 6)(4|5 + 6|1](6|4 + 5|3] ' 

where 

t = (K 1 + K 6 ) 2 . 

All other double-box coefficients can be computed by similar calculations. 

The calculation of the coefficient c 2 can be generalized for the configuration consid 
ered in fig. 18. 
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n-1 m n 

m+1 111 J 

Fig. 18: An infinite family of n-gluon double box scalar integrals. 
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There is only one helicity configuration that contributes to the octa-cut. A similar calcu- 
lation gives 



u st 



2 [1 2] [2 3] ... [to - 1 to] (m + 1 to + 2} ... (n - 1 n) 
1 



x 



(to + l\K m+1 , m+ 2,..., n \l}{n\K m+1 , m+2 ,...,n\m} 
where s and t are given by eq. ( |5.16| ) and u is given by 



U — (-Km+l,m+2,. 



{Km+1 + K m+2 + . . . + K r , 



(5.20) 



(5.21) 



6. Application to Three and Higher Loops 

Ideas presented in the previous sections can be applied to higher loops. Let us consider 
triple-box configurations appearing at three-loops. The configurations we consider are 
obtained from the double-box configurations at two loops by adding three new propagators 
to form the third loop. This way, one can produce a ladder diagram as in fig. 19a or a 
double box with a pentagon as in fig. 19b. We make a slight abuse of terminology and 
call both kind of configurations triple-box diagrams. 



(a) (b) 

Fig. 19: Three-loop triple box configurations, (a) A triple box ladder integral. 
(b) A double box with a pentagon. 



Every triple box contains ten propagatorsll. Therefore it is natural to start with a ten- 
particle cut. This produces ten delta- functions whereas the number of integration variables 

7 Of course, for a large enough number of gluons one can also find split triple boxes which can 
have 11 or 12 propagators. 
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is twelve. However, it follows from our previous analysis that box configurations develop 
additional propagator-like singularities which can also be cut (replaced by their discontinu- 
ities) . A triple-box configuration naturally admits two extra propagator- type singularities 
which allows us to consider twelve-particle cuts. Therefore, it should be possible to com- 
pletely localize all momentum integrals, at least if the number of gluons is not big enough. 
Then it is straightforward to write down an expression for the coefficients analogous to 



eq. ( 4.71) . Obviously, a similar analysis can be performed at any number of loops. It is 
interesting to mention that at three loops, some of the triple boxes that enter in the cal- 
culation of the four-gluon amplitude are not scalar triple boxes |[44|| . This means that the 
numerator in the integrand is not one but an inverse propagator. See fig. 20. This was 



also found to be the case for higher loops |44j] . 

We start our discussion with the analysis of ladder diagrams which allow a straight- 
forward generalization of our discussion in section 4. Then we turn to the triple-box 
integrals where one of the "boxes" has five propagators and see how our analysis of singu- 
larities realizes the phenomenon mentioned above. For concreteness, we will concentrate 
on the four-gluon amplitude though an identical analysis can be performed regardless of 
the number of external lines. Let us start with the triple-box configuration in fig. 21. 



a+hf 



Fig. 20: Schematic representation of a modified "triple box" integral given in [44] 
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K 2 



Fig. 21: A ladder triple-box configuration with four external gluons. 



+ 



r ■ _ q + _p+. 



- c + 



+ 



h/\ \/\ 



< 



+ 



< 



> 



1 \l\ 



+ 



r + q + p + 

2 n 2 *2 



Fig. 22: One of the helicity configurations contributing to the coefficient of the 
ladder triple-box in fig. 21. 



In this case, there are twelve helicity configurations, all giving the same answer. It is 
enough to consider one of them, for example the one in fig. 22. Similarly to the two- loop 
case, it is enough to consider the ten-particle cut because the measure integral factors out 
and cancels. Then the corresponding coefficient is given by 

4 = -i*irv? , {r l, T ,f 'V'r (6-1) 

(r 2 3) (3 ri) [4 r\[r r x \ 

where A^ 66 is the tree level four-gluon amplitude with the external lines (K± , K^-, r£, ft), 
s = (Ki + K 2 ) 2 and t = (Ki — r) 2 . In eq. (|6.1[), we used that the product of the six tree 



amplitudes was computed in section 3. Using momentum conservation, eq. (|6.1|) can be 
simplified to give 

di = -iAl ree s s t. (6.2) 



where t = (K2 + K3) . This coincides with the answer from [44 
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Fig. 23: A triple-box configuration. 



Now let us consider the configuration in fig. 23. Note that one of the loops has five 
propagator and this is why we said that the configuration was a double box with a pentagon. 
This is the basic reason why the integral is not a scalar box integral as we will see below. 
Let us start our analysis with the measure integral 

Jo = / dHd 4 qd 4 pS(l 2 )5((l - K 3 ) 2 )5((l -K 3 - K 4 ) 2 )5((l - K 3 - K 4 - K t ) 2 ) 

J (6-3) 

8((p -K 2 - l) 2 )5((p - K 2 ) 2 )5((q - Ki) 2 )5((p + q) 2 )S(p 2 )5(q 2 ) 

and perform the integration over p and q. After integrating over p we obtain (up to the 
momenta-independent factor) 

1 (6-4) 



(K 2 + q y{K 2 + iy 

Then we cut the "propagator" 1/(K 2 + q) 2 . This gives the fourth delta-function which 
allows us to perform the integration over q. This produces (again, we ignore the momenta- 
independent factors) one more factor of [K 2 + I) 2 in the denominator. Thus, we obtain 
that this triple box configuration has a singularity 

2 

(6.5) 



XK 2 + l)\ 

Let us calculate the coefficient of this triple box integral by calculating the product of the 
eight-gluon amplitudes. For every non-vanishing helicity configuration, the product of the 
six gluon amplitudes on the right is the coefficient c t of the two-loop four-gluon amplitude 
studied in section 3 and it is given by 

-A't ee st' 2 , (6.6) 
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where A'^ ree is the tree four-gluon amplitude with external lines (K ± , K 2 , / + , (/ — K3 — 
K±) + ) and t' is given by 

t' = (K 2 + l) 2 . (6.7) 

Therefore, in attempting to calculate the integral of the product of the tree amplitudes, 
the singularity (|6.5|) cancels out. This means that the coefficient of the triple box in fig. 23 



is zero. This is in agreement with results of [44|. On the other hand, the amplitude A'l ree 



has a factor 

(27) " (^W (6 ' 8) 

This indicates that the actual diagram has a singularity j^p^i. In order to account for 
this singularity we have to introduce a slightly modified triple box integral schematically 
shown in fig. 20. The basic idea is to multiply in the numerator by (K 2 + 1) 2 in order to 
cancel one of the power in (|6.5| ) and get the correct 1/(K 2 + I) 2 singular behavior. This 



shows that this triple-box integral should be in the list of scalar integrals of the amplitude 



under study. This is completely consistent with results of [44]. Now we can cut this 



"propagator" and completely localize the integral. Note that the combination 

[2l}S((K 2 + l) 2 ) (6.9) 

is not necessarily zero. This just means that we have to choose the solution 

Xi ~ A 2 . (6.10) 

The coefficient of this modified box is now straightforward to compute. The answer is 

d 2 = -iA\ ree s 2 t. (6.11) 

This coincides with the corresponding coefficient from MM. All other coefficients of this 



amplitude can be found in a similar manner. Even though we concentrated on the four- 
gluon amplitude, we could do the same analysis for any amplitude admitting additional 
cuts. 

As mentioned in the introduction, the basis of integral for M = 4 L-loop amplitudes 
of gluons is not known except for L = 1. One can imagine that a more systematic analysis 
along the lines of the discussion presented in this section might give a way of obtaining 
such a basis. It would be interesting to explore this direction in the future. 
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7. Conclusion 

In this paper, we observed that certain scalar double-box integrals which appear at 
two loops in M = 4 Yang-Mills theory possess hidden singularities. Such singularities are 
manifest after a quadruple cut is performed on one of the boxes. The end result is that 
one can straightforwardly calculate the coefficient of such integrals by an octa-cut which 
localizes the cut integral. The form of the coefficient is universal and it is given by the 
product of a certain number of tree-level amplitudes. This technique is applicable to all 
scalar double box integrals in amplitudes with less than seven external gluons and to a large 
subset of double box integrals for seven or more external gluons. The basis of integrals at 
two loops is not known in general. For four gluons the amplitude is given in terms of only 
scalar double-box integrals. If it turns out that the basis of integrals for five- and six-gluon 
amplitudes is also given by scalar double box integrals, then our technique gives a simple 
way of computing all those amplitudes for any helicity configuration. We also argued that 
this technique can be applied to higher loop amplitudes. At three loops we found that our 
technique can be easily extended to compute the coefficient of ladder diagrams. For the 
class of diagrams with a pentagon, our method shows that the coefficient of scalar integrals 
is zero and naturally gives the modified integral for which the coefficient does not vanish. 
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